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Abstract 

We have searched for molecular gas in the cD galaxy of a poor cluster of galaxies AWM7 using Nobeyama 
45 m telescope. We do not detect CO emission in the galaxy. Our limit of molecular gas in the inner 7.5 
kpc is Mh 2 < 4 x 10 8 M®. We estimate the total mass of molecular gas left in the cD galaxy when the gas 
deposited by a cooling flow once becomes molecular gas and the molecular gas is continuously evaporated 
by the ambient hot gas. The observational limit of molecular gas requires / > 10 -3 , where / is the ratio of 
the heat conduction rate to that of Spitzer. However, this contradicts recent X-ray observations showing 
/ < 10~ 5 . Thus, the non-detection of CO cannot be explained by the evaporation, and most of the cooled 
gas predicted by a cooling flow model may not change into molecular gas in the cD galaxy. Moreover, 
we estimate the evaporation time of molecular clouds brought to a cD galaxy through the capture of 
gas-rich galaxies and find that these clouds should not be evaporated if / < 10~ 3 — 10~ 4 . Therefore, the 
non-detection of CO in a cD galaxy could constrain the total mass of the molecular clouds brought into it. 

Keywords: Galaxies: clusters of — Galaxies: evolution — Galaxies: intergalactic medium — Galaxies: 
individual (AWM7) 
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1. Introduction 

The centers of galaxy clusters are usually dominated 
by very massive (~ 10 13 M®) galaxies. These galaxies 
are often called D or cD galaxies. The observations of 
cold gas ( < 10 5 K) give us the clues of the formation and 
evolution of the cD galaxies. 

Cold gas in cD galaxies has been investigated from the 
viewpoint of cooling flows. The cooling time of intraclus- 
ter medium (ICM) exceeds the Hubble time (~ 10 10 yr) 
in the most region of clusters (Sarazin 1986). However, 
around cD galaxies, the density of ICM increases and the 
cooling time decreases to t coo i < 10 9 yr. In the absence 
of heating, the gas is inferred to be cooling at a rate of 
M C f ~ 100M® yr" 1 (Fabian 1994). We will refer to M CF 
as a mass deposition rate from now on. Thus, the total 
mass accumulated around the cD galaxies would result 
in ~ 10 12 M® if the cooling occurred steadily at the rate 
over the Hubble time. Although many observers have 
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tried to detect the cooled gas mainly in massive cool- 
ing flow clusters A/cf > 100 M® yr -1 , most of them 
could not detect such a large amount of cold gas. Using 
recombination line luminosities, Heckman et al. (1989) 
estimate that the total mass of 10 4 K ionized hydrogen 
is less than 10 8 M®. Observations of the atomic hydrogen 
21 cm line limit the total mass of optically thin Hi to less 
than ~ 10 9 Mg, (Burns et al. 1981; Valentijn, Giovanelli 
1982; McNamara et al. 1990). CO observations limit 
the mass in clouds similar to Galactic molecular clouds 
to less than lO 9 " 10 ]^ (McNamara, Jaffe 1994; O'Dea et 
al. 1994; Braine, Dupraz 1994). Among the cD galaxies 
observed so far (~ 30), only one exception is NGC 1275, 
the cD galaxy in the Perseus cluster. The molecular gas 
of ~ io 9 ~ 10 M m has been detected (e.g. Lazareff et al. 
1989; Mirabel et al. 1989; Inoue et al. 1996), although it 
is smaller than the prediction of the cooling flow model. 
These observations may imply that most of the cooled 
gas becomes something other than molecular gas such as 
low mass stars, or the actual mass deposition rate, Mcf, 
may be reduced by some heating sources. 
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Before we move to investigate these possibilities, we 
should consider another scenario, that is, cooling flows 
actually exist but the molecular gas deposited by the 
cooling flows is continuously evaporated by the ambient 
hot ICM. The evaporation time of a molecular cloud is 
given by 

tcvap oc n c R c T l Qy i , (1) 

where n c and R c are the density and radius of a molecular 
cloud, respectively, and Ticm is the temperature of hot 
ICM (Cowie, McKee 1977). In relation (0), the satura- 
tion of heat flux is ignored although the following result 
does not change significantly even in the saturated case 
(see §4). White et al. (1997) investigate the data of Ein- 
stein Observatory and find the relations Ticm °c M^p 
and r coo i oc M C p 5 , where r coo i is the cooling radius. 
Thus, we obtain the relation: 

i ovap oc n c R 2 c M-°- 75 (2) 

On the other hand, if the molecular gas is accumulated by 
cooling flows and the age of cooling flows is much larger 
than t vap, the mass of molecular gas per unit volume is 
given by 

m mo \ oc (McF/^ ool )iova P oc n c R 2 c M^ r ° (3) 

Therefore, if n c and R c do not depend on Mcf too much, 
clusters with small Mcf should have large m mo \. In these 
clusters, we could find molecular gas. 

Since CO has been searched mainly in massive cool- 
ing flow clusters (Mcf > 100 M® yr^ 1 ), the observation 
of clusters with small Mcf is important because of the 
above reason. Moreover, we note here another impor- 
tance of searching CO in small Mcf clusters. Even if 
the actual mass deposition rate AfcF is reduced by some 
heating sources, molecular gas may be brought to cD 
galaxies. For example, the capture of gas-rich galaxies is 
another possible supply route of molecular gas into cD 
galaxies. In this case, the detection of the molecular gas 
would be easier for clusters with small Mqf- This is be- 
cause the X-ray emissions of these clusters are weaker 
than those of clusters with large Mcf , which means that 
the heating should be less effective in these clusters and 
that the cold gas would be less affected by the heating. 

We search CO in the cD galaxy NGC1129 at the center 
of a poor cluster AWM7 with relatively small mass depo- 
sition rate (M C f = 41M® yr _1 ; Peres et al. 1998). Note 
that AWM7 is one of the most closely studied clusters in 
X-ray. Ezawa et al. (1997) and Xu et al. (1997) find the 
metal abundance excess of the ICM at the center of the 
cluster. This suggests that the ICM is not well mixed and 
the cluster has not experienced violent cluster mergers 
at least recently. Thus, the molecular gas in NGC 1129, 
if exist, would keep intact. Moreover, the abundance ex- 
cess, especially in the cD galaxy (Xu et al. 1997), implies 
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that there has been star formation activity at the cluster 
center. The excess iron mass in the central region (r < 27 
kpc) is 8 x 10 8 Ms,. Assuming the 1 M® iron is ejected 
into ICM per 100 Ms of stars formed, the observation 
shows that ~ 10 11 M® of stars have been formed in the 
region. On the other hand, the present star formation 
rate of NGC 1129 is 0.04 M® yr" 1 within 1.57 kpc form 
the center (McNamara, O'Connell 1989). If the distribu- 
tion of stars in the galaxy is oc r -2 , the star formation 
rate for r < 27 kpc is 0.7 Ms, yr -1 . Thus, the present 
star formation rate in the region is smaller than that the 
average through the Hubble time (~ 10 Mg, yr _1 ), and 
the star formation in the past must be larger than that 
at present. If the 'starburst' occurred recently, molecular 
gas used for it would be left until present. In this paper, 
we assume H n = 50 km s" 1 Mpc" 1 throughout. 



2. Observations 

The 12 CO (J = 1-0) line was observed toward the cen- 
ter of NGC 1129 (a = 02 h 51 n T3. s 3, <5i 950 = +41°22 m 32 s ) 
with the 45-m telescope at Nobeyama Radio Observatory 
in 1999 March and May. The half-power beam width 
(HPBW) was 15", which corresponds to 7.5 kpc at the 
distance of NGC 1129 (z = 0.017325). The aperture and 
main beam efficiencies were ?7a = 0.40 and t/mb = 0.48, 
respectively. 

We used two SIS receivers that can observe two orthog- 
onal linear polarizations simultaneously. Martin-Puplett 
type SSB filters were used for image sideband rejection. 
The system noise temperature (SSB) including the atmo- 
spheric effect and the antenna ohmic loss was 400-600 K. 
As receiver backends 2048-channel wide-band acousto- 
optical spectrometers (AOS) were used. The frequency 
resolution and channel spacing are 250 kHz and 125 kHz, 
respectively. Total bandwidth is 250 MHz. Calibra- 
tion of the line intensity was made by the chopper- wheel 
method, yielding the antenna temperature (T£) corrected 
for both atmospheric and antenna ohmic losses. We used 
the main beam brightness temperature (T m b = Ta/^mb) 
in this paper. The telescope pointing was checked and 
corrected every hour by observing the 43GHz SiO maser 
emission in a late type star S-Per or W-And. The abso- 
lute pointing accuracy was better than 5" (peak value) 
throughout the observations. 



3. Results 

In figure 1, we present the CO (1-0) spectrum which has 
been binned by 20 km s _1 and has had baseline removed. 
The spectrum shows no significant CO (1-0) features ei- 
ther in emission or in absorption. The 3cr upper limit to 
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Fig. 1.. Spectrum of NGC 1129 in AWM7 (z 
velocity resolution is 20 km s — 1 



0.017239). The 



the cooling flow occurred steadily at the rate over the 
Hubble time. 

As mentioned in §1, when the molecular gas deposited 
by a cooling flow is continuously evaporated by the am- 
bient hot ICM, the detection of molecular gas would be 
relatively easy in clusters with small Mcf- Although we 
cannot detect CO, it constrains the evaporation rate of 
molecular gas and the heat conduction rate of the ICM. 
Using the results, we could investigate whether the non- 
detection of CO is consistent with the evaporation model. 

The accumulation time of molecular gas is 



iacc = M mo i/m< 1 x 10 8 yr 



(7) 



the flux integral is given by 
3<r ch AV 



'CO 



y/AV/AV ch 



K km s 



(4) 



where a c is the channel-to-channel rms noise, AT4h is the 
smoothed velocity channel spacing, and AV is the width 
of line. From figure 1, we obtain er c h = 0.006 K. 

We assume that the column density of molecular hy- 
drogen is 



7V H , = 2.8 x 10 20 /« 



CO cm 



(5) 



(O'Dea et al. 1994). The total mass of molecular hydro- 
gen is given by 



M mo i = 



4 In 2 



iV H2 m H2 



(6) 



where r is the beam size at the distance of the source 
and mHj is the mass of the hydrogen molecule (O'Dea 
et al. 1994). We assume that AV = 300 km s^ 1 for a 
rectangular line feature. This is the same as McNamara 
and Jaffe (1994) and nearly corresponds to the internal 
velocity dispersion of NGC 1129 (McElroy 1995). From 
equation (0), we obtain M mo \ < 4 x 10 s within 7.5 
kpc from the center. This is one of the most sensitive 
limits for cD galaxies. 

4. Discussion 

Although Mqf °f AWM7 is relatively small, the non- 
detection of molecular gas conflicts with a cooling flow 
model if the cooled gas becomes molecular gas and if we 
ignore the effect of the evaporation. Peres et al. (1998) 
estimate that the mass deposition rate and cooling ra- 
dius of AWM7 are M CF = 41M® yr -1 and r coo i = 
103 kpc, respectively. The analysis based on a cool- 
ing flow model shows that the mass deposited within r 
is M{< r) = M cv {r/r cool ) (Fabian 1994). Thus, the 
mass deposition rate within the beam of Nobeyama 45m 
is 77i — 3.8Mg, yr -1 , considering the projection effect. 
Thus, molecular gas of ~ 10 10 Af® would be detected if 



Although the properties of the molecular clouds de- 
posited by cooling flows are not well-known, we could 
calculate the evaporation time of the clouds as follows. 
In disk galaxies, molecular gas is considered to be pro- 
duced through disk instabilities (e.g. Larson 1987). In 
elliptical galaxies like cDs, we expect that the mechanism 
is ineffective. Instead, we expect that the molecular gas 
is produced through the thermal instability of ICM. One 
possible seed of the instability is supernova remnants (Fu- 
jita et al. 1996, 1997). Thus, we assume that a supernova 
remnant is the seed of a molecular cloud and that only 
thermal evaporation affects the cloud after the formation, 
although these may oversimplify the evolution of molec- 
ular clouds (refer to Loewenstein and Fabian [1990] for 
more realistic discussion about the issue). Note that the 
results in the following can be applied to other forma- 
tion mechanisms of cloud if the resultant mass is nearly 
the same. If we can ignore the fragmentation and coales- 
cence of molecular clouds, the mass of a molecular cloud 
is equal to that of a supernova remnant. Since the radius 
of a supernova remnant is given by 



R s ~ 50 pc 



Pi 



ICM 



4 x 10 5 cm~ 3 K 



-1/3 



10 51 ergs 



1/3 



(8) 



(Fujita et al. 1997), the mass of a molecular cloud is 



Mr. = 



4 „ 

-7ri? s 77i H n-iCM 



130 Mr, 



nicM 



P 



ICM 



4 x 10 5 cm -3 K 



10 51 ergs 



(9) 



(10) 



where Picm and tiicm are the pressure and the density 
of ICM, respectively, Esn is the energy released by a 
supernova, and mn is the hydrogen mass. 

If molecular clouds are in pressure equilibrium with the 
ambient ICM, the density of the molecular gas is 



n c = 4xl0 4 cm" 3 



P 



ICM 



4 x 10 5 cm~ 3 K 



10 K 



,(H) 
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where T c is the temperature of the molecular gas. Since 
wicmWh^ = n- c rri}i 2 R^, the radius of a molecular cloud 
is 



R c ~ 0.25 pc 



T 



Rs \ ( T c 
50 pc/ ^10K 

1/3 



1/3 



ICM 



4 x 10 7 K 



(12) 



We assume that heat is supplied to molecular clouds 
from the isothermal X-ray gas component prevailing even 
in the central region of clusters (e.g. Ikebe et al. 1999). 
Moreover, we assume that the isothermality of the com- 
ponent is retained by adiabatic heating or magnetic loops 
connected to the overall thermal reservoir of the cluster 
(e.g. Norman, Meiksin 1996) or the large heat conduc- 
tion rate of the gas constituting the component. The 
evaporation time of a molecular cloud embedded in ICM 
is given by 



1 x 10 5 yr 



1 



/ \A x 10 4 cm" 3 
xi^l l-^?-) (13) 



0.25 pc/ v 4x1 o 7 k 



1 



1 x 10° yr - 



/ \4 x 10 4 cm -3 



1/3 



Mc \ I 1 ICM 



130 Ma, J V4xl0 7 K 



(14) 



(Cowie, McKee 1977). Equation @ shows that if M c is 
larger than that given in equation (|l0|), t cvap should be 
larger. The parameter / is the ratio of the heat conduc- 
tion rate to that of Spitzer (1962) and < / < 1 (Cowie 
and McKee [1977] assume / = 1). When / < 1, the heat 
conduction rate and the mean free path of an electron are 
considered to be regulated by plasma instabilities around 
the cloud (Pistinner et al. 1996; Hattori, Umetsu 1999). 
If the mean free path of an electron is comparable or even 
greater than the radius of a cloud, the thermal evapora- 
tion is saturated. Defining the saturation parameter, 



a = 2700 / 



TlCM 

4 x 10 7 K 



V 10- 2 cm~ 3 ) 



Rc 



0.25 pc 



(15) 



the saturation occurs when do>l (Cowie and McKee 
1977). In the saturated case, the evaporation time is 
given by 



L cvap 



8 x 10 6 yr ' " 



x 



4 x 10 4 cm- 3 . 



Rc 



a0- 2 cm- 3 / V0.25pc 



ICM 



4 x 10 7 K 



-1/2 



2700 



8 x 10 6 yr 



™ICM 



4 x 10 4 cm- 3 



-3/8 



2/3 



(16) 



cm 



-3 



io- 2 

4 x 10 7 K 



1/3 



V130 M®, 

~ 1/2 ' do \ -'• * 



2700 



(17) 



(Cowie, McKee 1977). 

We will adopt T c = 10 K and E SN = 10 51 ergs; the 
temperature of a molecular cloud is the typical one in 
the Galaxy (Scoville, Sanders 1987). If we adopt the ob- 
served values nicM = 1-82 x 10- 2 cm- 3 and Ticm = 3.63 
keV (Mohr et al. 1999; Ezawa et al. 1997), equation 
(||) and (|l^) yield R c — 0.20 pc and equation ( |l5| ) shows 
that the saturation occurs when / > 5 x 10" . If the 
age of a cooling flow (to ~ 10 10 yr; Fabian 1994) is much 
larger than the evaporation time of a molecular cloud, the 
evaporation of molecular gas should be balanced with the 
accumulation. In this case, the evaporation time should 
be equal to the accumulation time. From equations (^|), 
(Gjl), (15), and (|l6|), this requires />1Q- 3 . However, 



recent observations of ASCA show that ICM is inhomo- 
geneous in temperature at least in some clusters (e.g. 
Ikebe et al. 1999). In order to explain this inhomogene- 
ity by the cooling flow model, / must be less than 10 -5 
at least around cooler X-ray gas (Pistinner et al. 1996; 
Hattori, Umetsu 1999) if the cooler gas component is not 
isolated by something like a magnetic field. Therefore, 
as long as / < IO -5 , the evaporation cannot account 
for the non-detection of CO and most of the cooled gas 
may become something other than molecular gas such as 
dust (Fabian et al. 1994; Voit, Donahue 1995; Edge et 
al. 1999) or low mass stars (Sarazin, O'Connell 1983), or 
there may be something wrong in the cooling flow model, 
that is, the actual mass deposition rate is much less than 
the one estimated by X-ray observations. 

So far, we have not considered the molecular gas 
brought by gas-rich galaxies merged into cD galaxies in- 
cluding NGC1129. Finally, we examine the evaporation 
of this kind of molecular gas. The mass of a molecular 
cloud brought through the capture would be larger in 
comparison with the case of cooling flows (~ 100 M®; 
equation p|). If molecular clouds in captured galaxies 
are similar to those in our Galaxy, the masses are typ- 
ically M c >10 5 M m (Binney, Tremaine 1987). Thus, if 
we adopt equation (|ll]) and the normalizations therein, 
i? c >2.5 pc. From equations ([if]) and (H), we expect 



t, 



cvap 



10 12 yr for / = 10" 



we obtain t 



cvap 



10 1 



yr- 



" 5 . Moreover, even if /= 10 3 , 
Thus, the evaporation can be 



ignored. This means that if molecular gas is brought into 
a cD galaxy through galaxy captures and /<10" 3 , the 
gas should be left there. (It is to be noted that when 
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Ticu ~ 10 keV, equation p3 shows that the condition 
£cvap ^ 10 10 requires / < 10 .) Hence, the non-detection 
of molecular gas strongly constrains the amount of molec- 
ular gas brought into the cD galaxy through galaxy cap- 
tures. Using a theoretical model based on a hierarchi- 
cal clustering scenario, Fujita et al. (1999) predict the 
amount of the molecular gas and compare it with the 
observations. 

5. Conclusions 

We have searched for CO emission from the cD galaxy 
NGC 1129 in AWM7. We have obtained the upper limit 
of molecular hydrogen mass (4 x 10 8 -M®). This is one 
of the most sensitive limits for cD galaxies. We predict 
the total mass of molecular gas left in the cD galaxy on 
the assumption that while the gas deposited by a cool- 
ing flow once becomes molecular gas, the molecular gas is 
continuously evaporated by the ambient hot gas. We find 
that the upper limit of molecular hydrogen mass shows 
/ > 10~ 3 , where / is the ratio of the heat conduction rate 
to that of Spitzer (1962). However, this is inconsistent 
with recent X-ray observations showing / < 10~ 5 . Thus, 
most of the cooled gas predicted by a cooling flow model 
does not seem to become molecular gas in the cD galaxy. 
Therefore, if / << 1 as is suggested by the X-ray obser- 
vations, the ultimate fate of most of the cooled gas may 
be something other than molecular gas such as dust or 
low mass stars. Alternatively, the actual mass deposi- 
tion rate may be much less than the one predicted by a 
cooling flow model. 

We find that molecular clouds brought to a cD galaxy 
by the gas-rich galaxies captured by the cD should not 
be evaporated when / < 10~ 3 — 10~ 4 . This implies that if 
we obtain the upper limit of the mass of molecular gas in 
a cD galaxy, we could constrain the supply of molecular 
gas brought into through the galaxy captures. 
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